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Abstract: Many dark matter models involving weakly interacting massive particles (WIMPs) 
feature new, relatively light pseudoscalars that mediate dark matter pair annihilation into 
Standard Model fermions. In particular, simple models of this type can explain the gamma 
ray excess originating in the Galactic Center as observed by the Fermi Large Area Telescope. 
In many cases the pseudoscalar’s branching ratio into WIMPs is suppressed, making these 
states challenging to detect at colliders through standard dark matter searches. Here, we 
study the prospects for observing these light mediator states at the LHC without exploit¬ 
ing missing energy techniques. While existing searches effectively probe pseudoscalars with 
masses between 5-14 GeV and above 90 GeV, the LHC reach can be extended to cover much 
of the interesting parameter space in the intermediate 20 - 80 GeV mass range in which the 
mediator can have appreciable Yukawa-like couplings to Standard Model fermions but would 
have escaped detection by LEP and other experiments. Models explaining the Galactic Cen¬ 
ter excess via a light pseudoscalar mediator can give rise to a promising signal in this regime 
through the associated production of the mediator with bottom quarks while satisfying ah 
other existing constraints. We perform an analysis of the backgrounds and trigger efficien¬ 
cies, detailing the cuts that can be used to extract the signal. A signihcant portion of the 
otherwise unconstrained parameter space of these models can be conclusively tested at the 
13 TeV LHC with 100 fb“^, and we encourage the ATLAS and CMS collaborations to extend 
their existing searches to this mass range. 
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1 Introduction 

Light, weakly-interacting massive particles (WIMPs) are a particularly compelling class of 
particle dark matter (DM) candidates. The case for WIMPs with masses close to the elec- 
troweak scale has been strengthened by recent observations of an excess in gamma rays origi¬ 
nating from the Galactic Center (GC) by the Fermi Large Area Telescope [1-10]. This signal 
has garnered much recent attention, since its morphology closely resembles that expected 
from dark matter pair annihilation into bottom quarks [9, 11], though other final states can 
also provide a good fit when systematics are properly taken into account [12]. Moreover, the 
signal suggests a WIMP annihilation rate close to that required in the early universe for a 
thermal relic to saturate the observed dark matter density [9], and the excess is difficult to 
explain in terms of astrophysical backgrounds alone [9, 13]. This has led many to believe that 
the Fermi GC signal may represent the first (indirect) observation of dark matter to date. 

A common and well-motivated class of models that can explain the observed excess 
features dark matter annihilating through a light pseudoscalar with Yukawa-like couplings 
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to Standard Model fermions [14-17]. For example, these states appear generically in two 
Higgs doublet models and their extensions [18], as well as pseudo Nambu Goldstone bosons 
associated with the spontaneous breaking of a new global symmetry [19-21]. Their couplings 
to Standard Model fermions can arise at tree- or loop-level (see e.g. Ref. [22] for an example 
with heavy vector-like fermions). Since they couple to the visible sector, such pseudoscalars 
can constitute a portal to the dark sector, mediating the annihilation of dark matter (DM) 
into SM final states [14, 23-28]. 

Understanding how dark matter interacts with the visible sector is a crucial part of the 
current dark matter program. Direct detection experiments [29-32] and the observation of a 
Standard Model-like 125 GeV Higgs with a small invisible decay width [33, 34] have severely 
constrained Z- and Higgs boson-mediated scenarios [21]. As a result, much recent work has 
been devoted to studying various possibilities for new mediator particles coupling weakly to 
the Standard Model degrees of freedom. Of these possibilities, pseudoscalars stand apart for 
several reasons. For one, they do not predict sizable spin-independent direct detection signals, 
in contrast with scalar and vector mediators. Furthermore, current collider constraints on 
new pseudoscalar particles are generally weaker than those on new scalar and vector degrees 
of freedom [35, 36]. 

If the GC excess is indeed a signal of dark matter annihilation, and if the annihilation is 
mediated by a new pseudoscalar particle, it is both important and timely to consider how one 
might probe such scenarios at colliders. Much progress has already been made on this front. 
Based on the topology and kinematics of the dominant dark matter annihilation channel, 
scenarios explaining the GC excess with pseudoscalar mediators can be grouped into roughly 
three types, each with distinct prospects for collider discovery: 

1. Models which rely on dark matter annihilating into on-shell mediators [37-42]. In 
this case, the annihilation rate into SM fermions factorizes and the coupling of the 
pseudoscalar mediator to SM degrees of freedom can be very small. Prospects for 
direct collider searches are often dim in this case, but there may be other handles on 
these models provided by direct detection, as well as fixed target and other precision 
experiments [37-40]. 

2. Scenarios featuring a pseudoscalar mediator with a significant invisible branching frac¬ 
tion [22, 25, 26, 43-46]. This results in distinctive missing energy signatures at the LHC 
which can be effectively probed by bb+MET^ mono-jet, and other existing and planned 
LHC searches, as studied in detail in e.g. Refs. [22, 25, 26, 43]. 

3. Scenarios in which the pseudoscalar mediator is expected to have a small branching 
fraction into dark matter particles [14-17, 27, 28]. This can occur when the coupling 
between the dark matter and the mediator is small relative to the coupling of the 
mediator to Standard Model degrees of freedom, or when on-shell decays of the mediator 
into WIMP pairs is not kinematically allowed. Such scenarios can be more difficult to 
probe directly at the LHC than case 2, since they lack a distinctive missing energy 
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signature [14]. In concrete models of this type, rare Higgs decays can be constraining, 
however the resulting limits can be straightforwardly avoided in many instances, as can 
limits from LEP, the Tevatron, and H-physics experiments (see e.g. Refs. [15, 16]). 
While a signal would arise in indirect detection experiments, it has been shown that 
the dark matter and mediator in this case might avoid detection elsewhere [14]. This 
rather grim scenario is appropriately known as “Coy dark matter”. 

In this study we will focus our attention on case 3 above, as it is a generic yet largely 
unconstrained possibility, as we discuss below. We will restrict our attention to light medi¬ 
ators, with masses below 90 GeV, as pseudoscalars with larger masses are already probed 
by existing LHC Higgs searches. Furthermore, light pseudoscalars are very attractive from 
the standpoint of the Galactic Center excess, since they can provide an efficient resonant 
annihilation channel for the light dark matter masses suggested by the signal and, in some 
cases, allow for a p-wave annihilation channel into pairs of mediators to drive down the relic 
abundance without violating constraints from dwarf spheroidal observations [21]. In this sit¬ 
uation, on-shell decays of the pseudoscalar to pairs of dark matter particles are suppressed 
and WIMP production at the LHC through the mediator will be negligible. Our strategy 
will be to extend LHC coverage to such scenarios by probing the light pseudoscalar directly 
through its interactions with the Standard Model degrees of freedom. The discovery of such 
a new particle would constitute a great step forward in our understanding of the dark sector 
and open up many possibilities for further study, including more dedicated experiments to 
probe its coupling to dark matter directly. 

As we discuss below, the GC excess can suggest an appreciable mediator coupling to 
down-type fermions. Consequently, we focus on the associated production of the mediator 
with a 5-jet or bb pair. We will assume that the mediator couples to Standard Model fermions 
with strength proportional to their mass, as in models with minimal flavor violation (MFV). 
We And that, for a signiflcant range of mediator masses and couplings consistent with the GC 
excess, a promising signal is predicted in the 1-2 b + a production modes, with a r^r~. We 
also explore the possibility of a ^ decays, which is more promising for low masses and 

likely features lower systematic uncertainties. Existing searches for pseudoscalars motivated 
by the Minimal Supersymmetric Standard Model (MSSM) and Next-to-MSSM (NMSSM) 
currently probe mediator masses down to 90 GeV and in the low-mass region between 5-14 
GeV. However, we And that coverage can be extended to pseudoscalars in the intermediate 
mass range (between 20-80 GeV), which are promising for explaining the GC excess and 
would have evaded detection by LEP. We encourage both ATLAS and CMS to expand their 
analysis to include this region. In this study, we detail the cuts and kinematic variables that 
can be used to reduce the large backgrounds and show the extent to which the parameter space 
in these models can be conclusively tested at the 13 TeV LHC with 100 fb“^ of integrated 
luminosity. We demonstrate this using a simplifled model and show the application of our 
results to the otherwise unconstrained parameter space of the NMSSM that is consistent with 
the excess (the NMSSM can be mapped directly onto our model). 


3 


It is important to emphasize that, although we will focus on pseudoscalars mediating 
dark matter annihilation consistent with the GC signal, our study can be applied much 
more generally to any model featuring light mediators with significant coupling to isospin- 
down Standard Model fermions. Since we assume that the invisible branching fraction of the 
pseudoscalar is small, our analysis of the predicted collider signatures does not depend on the 
pseudoscalar’s coupling to dark matter, nor on the nature of the dark matter itself. 

This study is organized as follows: in Sec. 2, we discuss the simplified model used for our 
analysis, its relationship to the GC excess, and the existing constraints on light pseudoscalars. 
The following section (Sec. 3), details the collider signatures of the new mediator, as well as 
the backgrounds and trigger efficiencies relevant for our analysis. Our results for the LHC 
discovery potential of light psuedoscalar mediators are presented and discussed in Sec. 4, 
with further details of the analysis contained in Appendices A, B, and C. We then apply 
these results to the NMSSM in Sec. 5, showing that the searches we propose here can cover 
much of the parameter space consistent with the excess and that is currently unconstrained 
by other experimental searches. Finally, we summarize and conclude in Sec. 6. 


2 A Simplified Model 


For our analysis, we follow Ref. [14] and consider a light pseudoscalar that couples to Dirac 
fermion dark matter, y, and to Standard Model fermions, with effective Lagrangian 


P ^ .^DM _ 5 
£int D X - i 




i=u,c,t ^ i=d,s,b, 


OdVi 


afafi, 


( 2 . 1 ) 


e,)L4,T 


where yi—mijv are the SM Yukawa couplings, with v — 174 GeV. We have assumed that the 
pseudoscalar a couples to the SM fermions with strength proportional to their masses. The 
pseudoscalar couplings to up- and down-type fermions are further assumed to depend on the 
overall scaling factors, Qu^d^ which we take to be the same for all down- or up-type fermions 
These factors appear e.g. in Two Higgs Doublet models (2HDMs) and their extensions; in a 
Type H 2HDM, — 1/^n = tan /3, where [3 is the ratio between the two SU{2) Higgs vacuum 
expectation values. With the addition of a singlet that mixes with the SU{2) doublets, the 
effective couplings become gu = cot (3 cos 9 and gd = tan/3 cos0, where 9 is the mixing angle 
between the SU{2) and singlet pseudoscalars. 

Note that Ref. [14] considered the case in which gd — gu = This situation is very diffi¬ 
cult to probe at colliders. Explaining the Fermi GC signal with gu — gd — ^ can require rather 
large values of unless the annihilation is quite close to the 5-channel resonance. Often 
in ultra-violet (UV) complete models, a sizable value for ^dm occurs together with low mass 
WIMPs in parametric regions featuring a large invisible branching fraction of the Standard 
Model-like Higgs [23], which is not observed. On the other hand, for pseudoscalar-WIMP 

^These assumptions need not be the case to explain the GC excess, and our results can be applied beyond 
this set-up by appropriately rescaling the pseudoscalar production cross-sections and branching ratios. 
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couplings that are not too large, the Galactic Center excess suggests enhanced couplings to 
down-type fermions, as we show below. This situation is much more promising from the 
standpoint of LHC searches and, in some cases, is not probed by existing searches. 


2.1 Explaining the Excess 

Given the Lagrangian in Eq. 2.1, the zero-temperature s-channel annihilation rate for dark 
matter through a pseudoscalar into SM fermion pair fifi is 



( 2 . 2 ) 


where m^, Nc^i are the mass and color factor of the decay states, gi is either gu^d depending 
on the fermion, and is the total width of the mediator. Throughout this study we will 
assume that the dominant DM annihilation channel is XX tfb. This mode has received the 
most attention in explaining the GC excess, and although a recent analysis has pointed out 
that other channels can also explain the signal [12], annihilation into a bb pair still provides 
a very good fit to the data. 

There have been several recent developments in determining which annihilation channels, 
WIMP masses, and annihilation rates best fit the Fermi data. For the bb channel, most 
previous analyses had suggested that should fall roughly in the range 35 GeV < < 

50 GeV with annihilation rate {av) ^ 2-6 x 10“^^ cm^/s [9, 14] (the required annihilation 
rate for self-conjugate dark matter would be reduced by a factor of two relative to these 
values). However, there are large systematic uncertainties associated with the propagation 
of gamma rays in the Galactic Center that must be taken into account. The impact of these 
systematics was first studied in Ref. [13], and subsequently by Ref. [12], which performed a 
detailed analysis incorporating several different models for the diffuse gamma ray background 
supplied by the Fermi collaboration. The end result is that the range of WIMP masses and 
annihilation modes statistically consistent with the excess increased significantly once these 
systematic uncertainties were taken into account. In particular, the range for WIMP masses 
annihilating primarily into bb was extended to [12] 


(2.3) 


35 GeV < < 165 GeV, xx ^ bb. 


for specific values of the annihilation rate. 

Across this mass window, the signal from the Galactic Center suggests a clear range of 


values for the coupling constants of the mediator to SM states for a given ^dm in this setup. 


Since the annihilation cross-section and pseudoscalar width are dominated by down-type 
interactions (for BR{a xx) 1)? the only significant parametric dependence is on 
gd^ and m^. The down-type scale factor gd required to explain the GC excess for = 45 
GeV and = 145 GeV (close to the best fit mass for the Fermi model (d) from Ref. [12]) is 
shown by the bands on the left and right hand side of Fig. 1, respectively, as a function of 





rria [GeV] 



20 40 60 80 100 120 


Ma [GeV] 


Figure 1. Value of the mediator coupling to down-type SM fermions (relative to that of a Standard 
Model-like Higgs boson of the same mass) required to explain the GC excess as a function of rua, for 
several values of Qdm and masses = 45, 145 GeV, assuming that y saturates the relic abundance. 
The shaded regions are compatible with the signal, with the red (upper) regions in each band excluded 
by the recent dwarf spheroidal constraints from Fermi [48], and the yellow (lower) regions with an 
annihilation rate compatible with both the excess and the constraints. The upper bound on from 
existing LHG searches for a is shown in blue. 


for various values of ^dm- The range of annihilation rates allowed in the low mass case (LHS) 
is taken from Ref. [13], while the allowed values in the high mass case (RHS) are taken from 
Ref. [12]. In both cases the local dark matter density is assumed to be pQ = 0.4 GeV/cm^. 
The preferred regions depend on = J/Jo, the ratio of the angularly-averaged integral over 
the line-of-sight dark matter density Pdm (^)5 given by 

^(V’) = ^/' f dspDuirf, (2.4) 

J zXO J 1.0.s. 

to the canonical value Jq. For the low mass (m^ = 45 GeV) case we take J — while 
for — 145 GeV we take J — 0.3, which is within the systematic uncertainties discussed 

in Ref. [12]. The latter choice allows for an annihilation rate close to the canonical thermal 
freeze-out value {{(Jv) ~ 4.4 x 10“^^ cm^/s for Dirac fermion dark matter [49]) and consistent 
with the Fermi signal while evading the dwarf spheroidal constraints, discussed below. 

For reasonable choices for the value of must be quite large to account for the GC 
excess, unless the masses are tuned to fall very close to the resonance. In addition, reducing 
the X abundance has the effect of increasing the preferred value of g^ for a given ^dm- The 
regions of parameter space with large in many cases preferred by the signal, predict a 
significant mediator production cross section at the LHC in association with bottom quarks. 
Also, the pseudoscalar’s invisible branching fraction is small across the entire parameter space, 
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except for low and large ^dm- For '^a < 2m^ an on-shell pseudoscalar cannot decay to 
a pair of WIMPs, while for rria > 2m^ we find that BR{a xx) >0.1 only for ga ^ ^ for 
^DM = 0.1 in the = 45 GeV case, since everywhere else gd (^dm) is too large (small) for 
this decay to contribute appreciably to the total width. 

It is important to note that the Fermi collaboration recently released updated limits on 
the dark matter annihilation rate from observations of dwarf spheroidal (dSph) galaxies [48]. 
The resulting constraints^ are in mild tension with a dark matter explanation of the excess, 
however there is still a large amount of parameter space capable of explaining the GC excess 
that survives this constraint. This is shown in Fig. 1, in which the red bands show the impact 
of the dwarf spheroidal limits (points in these bands could potentially explain the excess but 
are excluded at 95% C.L.). Meanwhile the yellow bands show points consistent with both the 
GC excess and dSph limits. Note that in the high mass case all points consistent with the 
excess are compatible with the dSph constraints for our particular choice of J. 

One concern may be that, since the recent dSph constraints disfavor larger annihilation 
rates, some points with light WIMP masses consistent with the GC excess and dSph limits 
will tend to produce too large a relic abundance. The dark matter relic density is set by the 
annihilation rate at finite temperature, which can differ from that at T = 0. In particular, 
for s-channel annihilation through a pseudoscalar with rria < 2m^, the annihilation rate at 
T = 0 is greater than that at freeze-out {Tf^o^ ^ m^/20). The upper limit on the annihilation 
rate, set by the Fermi dSph results, is below the required annihilation rate at freeze-out 
for < 100 GeV, naively disfavoring this region. However, there are several well-known 
and straightforward exceptions to this reasoning [53]. For example, p-wave processes with 
contributions to the total annihilation rate scaling as (with vdm the relative dark matter 
velocity) will become important at freeze-out, increasing the annihilation rate at Tf^o. but not 
altering the T = 0 prediction. An example of such a process generically expected along 
with light mediators is xX (this is another virtue of the light pseudoscalar scenario). 

Other scenarios allowing for an enhanced annihilation rate at Tj.o. relative to that at late 
times include those with additional co-annihilation channels or featuring rria > 2m^ so that 
{(7v)t=o < {(Tv)T=Tf o • Thus, although in some cases the dSph limits may result in requiring 
some additional tuning or model-building to achieve the correct DM relic abundance, dark 
matter explanations of the excess, particularly those involving s-channel annihilation through 
a relatively light pseudoscalar, are alive and well. Note that this tension largely disappears 
above ^ 100 GeV, since the dSph upper bound is above the canonical WIMP cross- 
section in this region (although one should verify that contributions to the annihilation rate 
at freeze-out from the other states in the theory do not over-dilute the relic density). 

In summary, dark matter annihilating through a relatively light pseudoscalar can explain 

^ There are also potential constraints on dark matter interpretations of the excess coming from observations 
of the cosmic ray antiproton spectrum [50]. These constraints depend sensitively on the propagation of the 
charged cosmic rays through the galaxy, which is difficult to model and results in large uncertainties on the 
predicted flux [51]. With conservative choices for the propagation model, the signal can be shown to be 
consistent with the current limits [50] . Similarly, radio observations of the Galactic Center region can also be 
consistent with a dark matter interpretation of the excess [52] . 
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the Galactic Center excess and be compatible with the recent dwarf spheroidal limits from 
Fermi. In all most discussed and shown above we expect BR{a xx) 1? either because 
rria < 2m^, qdm 1, or both. This implies a low likelihood of observing the pseudoscalar 
through missing energy signals at the LHC. In the following subsection, we describe some 
of the other existing constraints on the parameter space and highlight the need for direct 
coverage of these scenarios at the LHC. 


2.2 Existing Constraints 

Our goal will be to ascertain to what extent LHC searches can cover the parameter space 
shown in Fig. 1 that is not currently probed by LHC searches [54-58] . To our knowledge, 
there are currently no direct constraints on the parameter space of our simplified model with 
15 GeV < rria ^ 90 GeV. By this, we mean that there do not exist constraints depending only 
on the pseudoscalar’s coupling to SM fermions in this mass range. There are indeed several 
other indirect constraints, but these are inherently dependent on other degrees of freedom in 
the UV complete theory and can be straightforwardly avoided in many cases. We will present 
explicit examples of points evading all of the searches discussed below but still predicting an 
observable LHC signal in the NMSSM in Sec. 5. 

Pseudoscalar mediators with GeV-scale masses predict highly suppressed direct dark mat¬ 
ter detection cross-sections. At tree-level, the pseudoscalar only interacts spin-dependently 
with nuclei. Using the expressions and results found in Refs. [14, 59], we find that the spin- 
dependent scattering cross-section for dark matter off of nuclei via the pseudoscalar is far 
below the reach of current and planned experiments (ctsd ^ 10“^^ cm^) across the parameter 
space we consider. This is thanks to the 1/m^ suppression in cjsd in this regime. Also, while 
spin-independent scattering can occur via one loop diagrams, this contribution is also much 
too small to be observed. The difficulty in observing dark matter interacting with the visible 
sector primarily through a pseudoscalar in direct detection experiments in indeed one of the 
main reasons such models are understood to be coy. 

Light pseudoscalars can also be constrained by flavor observables^. Loop diagrams in¬ 
volving the pseudoscalar can generate effective flavor-changing vertices [21, 61]. The limits are 
severe for pseudoscalars lighter than the B and T meson scale simply because the mediator 
can be produced on-shell in decays. For rua ^ 10 GeV, the constraints are very significantly 
relaxed, with the most stringent arising from LHCb [62] and CMS [63] measurements of 
BR{Bs For ^ tub ^ the limits are approximately 


9d ^ 


3ma 
10 GeV 


(2.5) 


^Other precision tests, such as the pseudoscalar contributions to the muon g — 2, do not impact the 
parameter space we consider [60], although a light pseudoscalar with very large gd can, in some cases, help 
reconcile the observed {g — 2)^ with the SM prediction [28]. Also, precision electroweak measurements are 
typically unconstraining, since pseudoscalar contributions to the gauge boson self-energies first appear at two 
loops. 




considering the pseudoscalar contribution alone [21]. This constraint would naively appear to 
directly constrain some of the parameter space shown in Fig. 1, however, the new contributions 
to Bs iJL~ are strongly model-dependent [64]. For example, in supersymmetric UV 

completions of our model, such as the NMSSM, there are several new contributions which 
enter with opposite sign to that from the a-induced vertex. Thus, cancellations can occur over 
large portions of the parameter space allowing for light pseudoscalars with large couplings to 
SM fermions (i.e. above the naive upper bound of Eq. 2.5) [16], once again highlighting the 
need for direct probes of this parameter space. 

For light mediators with 2ma < rrih {h is the 125 GeV SM-like Higgs), exotic Higgs 
decays to pseudoscalar pairs can affect the Higgs width and signal rates [65, 66], which are 
constrained by both ATLAS [67] and CMS [68]. Evidence for h ^ aa decays was also searched 
for at LEP [69] and the Tevatron [70]. Such decay modes can also be very effectively probed 
at the High Luminosity LHC [66]. Indeed, this has long been recognized as an important 
potential discovery channel of NMSSM pseudoscalars at colliders [71-75]. However, these 
constraints depend on the haa coupling which, in some cases, can be made appropriately small 
in realistic models [16, 23], especially those in which the pseudoscalar coupling to Standard 
Model fermions does not arise through mixing with the SM-like Higgs [22]. Alternatively, 
simply taking rria > rnh/2 avoids these constraints altogether. 

Another indirect constraint arises from LEP searches for e^e“ ^ ha production [76]. 
While these results prohibit MSSM-like pseudoscalars lighter than 90 GeV for all values of 
tan/3, these bounds depend on the Zha coupling, which is model-dependent, and can again 
be straightforwardly avoided [16]. For example, in a Type H 2HDM with an additional singlet 
(2HDM-hS), the Zhia coupling scales as 

gzhia ^ cos e {Sii sin /3 - Si 2 cos /3) (2.6) 

where Sn and Si 2 are the corresponding entries in the matrix diagonalizing the 3x3 CP- 
even mass matrix with the Higgs bosons ordered in mass (see e.g. Eq. 2.22 in Ref. [77]). 
Contrasting to ^ cos0tan/3, we see that the simple limit cos0 1, tan/3 ^ 1 can result 
in an appreciable with a significantly suppressed gzha- 

Finally, existing MSSM Higgs boson searches at the Tevatron [78-82] and the LHC [54, 55, 
57] constrain g^ for > 90 GeV, but in an effort to avoid the large backgrounds encountered 
for lighter masses, and because LEP had already ruled out MSSM-like pseudoscalars with 
masses below 90 GeV, the published limits do not extend below the Z mass. There are also 
searches for light (m^ < 15 GeV) pseudoscalars at CMS [56], motivated by certain limits of 
the NMSSM. However, the 15 GeV < rria ^90 GeV mass range remains currently untested 
Although the collider limits on a light pseudoscalar can be avoided, one might also be 
concerned about the consistency of this scenario once the model is UV-completed. Our 

^There are also searches for a ^ 77 that probe masses down to iria = 65 GeV [58], however, as we will see 
below, the production cross-section in the diphoton mode is very small in the parameter space we consider 
and thus significantly below the existing limits. 
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Lagrangian is not invariant under SU{2)l x C/(l)y, and so, given a particular UV completion, 
one should also check that constraints on the other states can be satisfied while demanding a 
light pseudoscalar. In 2HD+S models, for example, most constraints on the rest of the Higgs 
sector can be satisfied by simply taking the charged Higgs mass to be moderately heavy 
(a few hundred GeV) with an appropriate choice of tan/3 [54, 55]. Such requirements are 
consistent with light pseudoscalars and sizable as shown in e.g. Ref. [16] and in Sec. 5 for 
the NMSSM. 

Perhaps surprisingly, there is a significant gap in coverage for light pseudoscalars with 
appreciable couplings to SM fermions, as arise in models explaining the GC excess or oth¬ 
erwise. This situation has room for improvement. In the remaining portion of this paper, 
we will investigate to what extent searches similar to those already existing for heavy MSSM 
Higgs bosons and for light NMSSM pseudoscalars can directly probe the parameter space 
motivated by the Galactic Center excess. This task requires a careful treatment of the back¬ 
grounds below the Z mass. As we will show below, the backgrounds can be substantially 
reduced by using a suitably chosen sequence of kinematic cuts. 

3 Light Mediators at the LHC 
3.1 Production and Signals 

Heavy neutral Higgs bosons in two Higgs doublet models are being searched for via a variety of 
experimental signatures, including gluon fusion (ggF) production, or production in association 
with top or bottom quarks [54, 55]. These canonical Higgs-type searches become much more 
difficult below the Z threshold, where the backgrounds increase dramatically. Fortunately, 
as we have shown in Sec. 2 above, light pseudoscalar mediators consistent with the Galactic 
Center excess can have enhanced couplings to down-type Standard Model fermions relative 
to those expected for a Standard Model-like Higgs boson of the same mass. This results in 
an enhanced production cross-section in modes involving b quarks, and (potentially) in the 
gluon fusion channel relative to the Standard Model-like case. This situation is depicted on 
the left hand side of Fig. 2, which shows as an example the enhancement of both the inclusive 
bba (black) and gluon fusion (red) production cross-sections with = g~^ (i.e. cos0 — 1 ), 
relative to those with g^ = gu — ^ (^ 0)5 as a function of g^^. The enhancement of the 
bba cross-section is independent of m^, as it only depends on g^ for a given m^, while the 
differently styled red curves correspond to crggF/crggF,o foi* different values of the pseudoscalar 
mass. The enhancement is substantially larger in the bba mode across the parameter space, 
which suggests focusing on production processes involving b quarks rather than the gluon 
fusion process. 

We consider the branching ratio of the pseudoscalar into various final states, assuming 
BR{a xx) is negligible, on the right hand side of Fig. 2. The pseudoscalar’s branching 

^This relationship between Qu and ga holds in the MSSM, but not the NMSSM. In the latter case, Qu = 
cos^ which results in values of Qu that suppress the top quark contribution to the ggF loop process even 
further than we have considered here. 


10 



cos 0=1 



Figure 2. Left: The enhancement of the inclusive bba and gluon fusion production cross-section 
relative to those with = g~^ = 1 as a function of gd- The dotted, dash-dotted, dashed, and 
solid red lines correspond to the enhancement in ggF production for rria = 20, 40, 60, and 80 GeV 
respectively. The corresponding enhancement for bb associated production is shown by the solid black 
curve (the enhancement is independent of rria). Right: Branching fraction of the pseudoscalar into 
various final states (assuming BR{a xx) is negligible). Note that the branching ratios into fermions 
are nearly independent of gd (since the total width is set primarily by a ^ 66, decays), while 

the a ^ 77 partial width is substantially suppressed for gd> 


fraction into photons is small and is further suppressed for > 1 which, when combined with 
the increased backgrounds for rria < mz^ suggests that diphoton searches will likely be unable 
to probe the low-mass pseudoscalar mediators we are interested in. On the other hand, while 
the favored decay is into a 66 pair, searches for such resonances would contend with large, pure 
QCD backgrounds to exploit this mode. Thus, to avoid large backgrounds while maintaining 
a reasonable signal, and to maximize the enhancement of the production cross-section, we 
propose a search for the pseudoscalar in second and third generation dilepton and 

/i+/i“) pair production in association with one or two 6-jets. Of course this strategy requires 
that the pseudoscalar couples to leptons, which is typical in extended two Higgs doublet 
models, but need not be the case [21]. 

Similar searches have been considered by both ATLAS [54] and CMS [55] , but are focused 
on higher mass resonances motivated by two Higgs doublet models and the MSSM, where 
the mass region of interest is greater than about 90 GeV [76] due to LEP searches and 
precision constraints on heavy Higgs bosons. Also, previous theoretical studies in the context 
of the NMSSM have investigated the potential for the LHC to probe light pseudoscalars with 
somewhat similar searches [83-88]. However, these investigations did not incorporate trigger 
and detector effects, and did not analyze the effects of cuts on the signal and backgrounds 
in detail, which is a major component of this work and crucial for obtaining an observable 
signal. While Ref. [87] arrives at largely negative conclusions regarding bba production (at 
least in the NMSSM with partial universality), our analysis suggests a much more positive 
picture once appropriate cuts are implemented. 

It is worthwhile to point out that the CMS search in Ref. [56] hnds sensitivity down to 
gd ^ 3 for masses up to rria ^ 14 GeV in the gluon fusion mode with a One might 
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be inclined to conclude that this search channel could simply be extrapolated to larger masses 
in the scenarios of interest. However, this is unlikely to be the case. Fig. 2 shows that the 
gluon fusion production cross-section is actually suppressed for 1 < < 10 as compared to its 

value with = 1, given our assumptions about the couplings. The suppression increases with 
rria and is due to the decreased top quark loop contribution that is otherwise dominant for 
heavier masses. In addition, due to the kinematic beta factor y^l — ( 7 ^)^, the bb branching 
ratio is suppressed for smaller values of rria^ resulting in an increase in the ii~ branching 
fraction. For example, BR{a is enhanced by almost a factor of 2 at rria — 10 GeV 

versus rria > 20 GeV. Thus, for the scenarios we consider, production modes involving down- 
type fermions at tree-level would appear more promising than those relying on gluon fusion 
production and decays to muons, although different assumptions about the coupling structure 
could alter this conclusion. For a related analysis of the potential LHC reach in the Ob mode 
in Z' models, see e.g. Ref. [89]. 

In the remainder of this section, we discuss the challenges and strategies for examining 
low mass pseudoscalars with enhanced couplings to down-type fermions, > 1. 

We implemented our simplified model in FeynRules 2.0 [90], and generated both our 
signals and backgrounds at leading order (LO) using MadGraphB+aMCONLO [91]. We then used 
Pythia 6.4 [92] to decay the r leptons and hadronize the 6 -jets, and incorporated initial and 
final state radiation, with an appropriate scale used for the MLM matching of hard element 
and radiated jets. Detector simulation for trigger and tagging was performed using Delphes 
3.0 [93]. Trigger effects were implemented as step-function cuts at the analysis stage, though 
some minimum kinematic requirements were enforced at the generation phase. 

Diagrams for some of the primary production modes for the signal are shown in Figure 
3. To avoid the appearance of potentially large logarithms arising from the phase space 
integration over collinear final state quarks, the semi-inclusive b(b)a events were generated 
with b quarks included in the parton distribution functions (pdfs) of the proton. This is known 
as the “five fiavor scheme” which effectively re-sums the large logarithms [94-96]. Exclusive 
bba events were generated without the inclusion of the b pdfs since the resulting contributions 
are doubly pdf-suppressed and subleading when compared to the gluon induced processes. 
To avoid double counting between the two-body, 6 ( 6 )a, production and the three-body, 66 a, 
production mode where one of the 6 -jets is collinear with the proton beam, the three-body 
production mode was generated with a minimum > 5 GeV. 

There are several technical difficulties associated with accurately calculating the two- 
body 6 ( 6 )a production cross-section at hadron colliders, which have received much attention 
in the literature [97-104]. In particular, the leading order production cross-sections are known 
to exhibit a substantial dependence on the renormalization and factorization scales, pr smd 
/ij, respectively [99, 103]. For our signal generation, we consider dynamic scales defined by 

Mr = M/= { E (3-1) 
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Figure 3. Some of the diagrams contributing to the production of the pseudoscalar, a, at the LHC. 
The two rightmost diagrams arise in the 5FS. 

where / is an overall scaling factor, and i refers to the produced 6’s and a. This is in keeping 
with previous analyses in the context of Standard Model-like Higgs production [99, 104-107]. 
We considered the impact of the scale dependence by varying the overall scaling factor in the 
range [1/2,2], which resulted in a 2-20% change in the production cross section, with larger 
effects occurring for smaller values of m^. This is consistent with the range typically found 
in the literature [98, 99, 101]. 

To further validate the results of our leading order calculation, we have compared our 
LO result for the dominant {gb{b) b(b)a) production mode to the next-to-leading order 
(NLO) result calculated in the five flavor scheme implemented in the program MCFM [108] for 
several choices for (we neglect the difference between scalar and pseudoscalar production 
which are small [104]). We find that our LO results exhibit reasonable agreement with the 
NLO result, falling within a factor of 1-2 across the parameter space we consider. Addition¬ 
ally, there are theoretical uncertainties related to the specific choice of parton distribution 
functions, which have been shown to be of order ^ 5% for low masses [103], as well as some 
residual renormalization scheme dependence (MadGraph uses an on-shell scheme, while e.g. 
MCFM uses MS). To account for these effects, Appendix C takes a conservative approach and 
explores the effect of a factor of 2 over-estimation in our signal and, separately, a factor of 
2 under-estimation in the backgrounds. Our overall conclusions are not significantly affected 
by this re-scaling, and so we believe them to be quite robust. 

For an experimental search, we consider three possible leptonic tagging channels: SRI 
requires one electron and one muon; SR2 requires one lepton (e or jj) and one hadronic r; SR3 
requires two muons. SRI is motivated by excellent trigger response, while SR2 is motivated 
by the larger branching ratios and SR3 is motivated by a resonance search methodology in 
the di-muon invariant mass spectrum that allows for the use of data-driven backgrounds. In 
all three signal regions, we also require 1-2 6-jet tags, and no light jets, where light jets are 
defined as pT > 40 GeV. The signals are therefore inclusive for light jets with pT < 40 GeV, 
such as those that are commonly generated from ISR effects. These tagging requirements 
significantly suppress fake backgrounds arising from vector boson production in association 
with light jets. 
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We assume the default CMS tagging efficiencies that are implemented in Delphes 3.0, 
which are as follows. For tagging, electrons are required to have pT > 10 GeV and \r]\ < 2.5. 
Within the inner region of the detector, \r]\ < 1.5, we assume a tagging efficiency of Cg = 0.95, 
while for the outer region but with \r]\ < 2.5 we assume Cg = 0.85. The rate at which jets 
fake electrons is taken to be /g = 0.0001 and uniform over the whole detector. For muons, 
we require that candidates have pT > 10 GeV and \r]\ < 2.4. Since our analysis involves 
only low pt muons, we take a fixed tagging efficiency of = 0.95, which is appropriate for 
pi^ < 1000 GeV. For the tagging of hadronic taus, we require I 77 I < 2.5 and take a fixed tagging 
efficiency of e^- = 0.4 with a fake rate for mistagging a light jet as a hadronic tau of ^ = 0.001. 
The tagging of 6 -jets occurs only where p^ > 15 GeV and \r]\ < 2.5, with an efficiency of 
= 0.5 tanh(0.03pr —0.4) within the inner detector, \r]\ < 1.2, and = 0.4 tanh(0.03pr —0.4) 
up to the boundary of \p\ < 2.5. Light jets are taken to fake 6 -jets at a rate of ^ = 0.001, while 
c-jets faking 6 -jets follow a formula similar to the 6 -tagging efficiency but with coefficients of 
0.2 and 0.1 for the two regions, respectively. 

3.2 Trigger Effects 

Since the signal typically produces very soft jets and leptons, trigger effects are very important 
to consider. To account for the effect that trigger has on our results, we have implemented 
a variety of triggers as a step-function cut based on what we believe are reasonable off-line 
triggers for CMS^. The following primary triggers are potentially relevant to our study: 

• le: single electron with pT > 35 GeV; 

• 1 /i: single muon with pT > 25 GeV; 

• 2 /i: di-muon leading with pT > 17 GeV, subleading pT > 10 GeV; 

• CTh'. electron -h hadronic tau with pj. > 45 GeV, p^ > 19 GeV; 

• jiTh'- muon -h hadronic tau with pTj. > 40 GeV, p^ > 15 GeV; 

• e/i: leading electron -h muon with p^ > 23 GeV, pi^ > 10 GeV; 

• pe: electron -h leading muon with p^ > 12 GeV, p^ > 23 GeV; 

We also include other triggers, such as those involving photons, jets, plus MET, and 6 -jets, 
but these provide a negligible effect on the signal events (i.e. < 0.3% of signal events pass 
all the non-primary triggers combined) and so are not included in the above list. The non¬ 
primary triggers pass a significant portion of the backgrounds, however, which necessitates 
their inclusion, but this indicates that these events have distinctive signatures that can be 
eliminated from the analysis by kinematic cuts. 

^Dilepton triggers are motivated from discussions with James Hirschauer of the CMS collaboration. Final 
trigger details for LHC13/14 are not currently available. 
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Due to the low mass of the pseudoscalar in our search, a significant number of the 
production events will not pass the trigger. Since we are not privy to the details of the final 
triggers, we consider the effect of varying the muon pT thresholds for the triggers that include 
a primary muon. These triggers have the greatest likelihood for discretionary variation in a 
dedicated experimental search, and are the most important due to having the lowest inclusive 
cross sections and thus pr thresholds. We analyzed the cross section of signal events that 
pass each of the primary trigger cuts ^ fraction of the cross section of generated 

events {cr^^Rx ~ x BR{t^t~ SRx) x esRx) with the same tagging signature in each 
signal region, independently: 


T 

rpTy _ ^SRx 

^SRx ~ ^gen 
^SRx 


(3.2) 


where SRx refers to the signal region and Ty refers to the specific trigger. This ratio can 
be considered as a sort of trigger efficiency. Of note, we found that the e + Th and fi -\- 
triggers did not pass any of a preliminary 200k generated events, likely due to the hard cut 
on the pt of the and the low mass of the pseudoscalar. Since the hadronic tau has a large 
fake background from mistagged light jets, we do not anticipate that the trigger threshold for 
hadronic tau pt will be improved enough to make these triggers worthwhile to consider. While 
the fake rate of jets for electrons is smaller than for hadronic taus, we believe it is unlikely 
that any significant improvement in the electron trigger thresholds will be implemented as 
there would still be a larger increase in the inclusive cross section than for similar changes in 
the muon trigger thresholds. 

The summary of the trigger efficiency ratios in Eq. 3.2 for the default implemented 
triggers is shown in Table 1, while an analysis of the effect of varying the threshold for the 
muon Pt in the 1/i, 2/i and pe triggers for each of the three signal regions is shown in Figure 
4. A naive interpretation of this figure suggests that the single muon trigger includes a larger 
fraction of the signal than pe or 2p triggers, but it is important to note that the single muon 
inclusive cross section at the LHC is significantly larger than the muon+electron or dimuon 
inclusive cross sections, and thus will typically have a higher pT threshold than the other 
triggers and a lower trigger efficiency, as shown in Table 1. 


3.3 Backgrounds and Their Reduction 

Since the QCD backgrounds at the LHC are significant, the fake rate of jets as electrons, 
hadronic r-jets, and 6-jets are important to take into account. Additionally, backgrounds with 
similar kinematics to the signal we examine produce soft leptons that may not be identified as 
easily or may fall outside of the central region of the detectors where tagging is possible. Thus, 
backgrounds producing more than two leptons, where one is not tagged, may contribute to the 
signal regions. To account for these effects, we include backgrounds that produce between one 
and three leptons (e, /i, r), and 0-2 6-jets, in association with 1-3 light jets (with nij + rij < 3), 
since our signal is inclusive to low pT light jets. The following background processes are 
generated: 
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SRx 

rua (GeV) 

le 

Ip 

2/i 

e + p 

p + e 

all 


20 

0.03 

0.12 

0.00 

0.16 

0.13 

0.30 

SRI 

40 

0.05 

0.21 

0.00 

0.22 

0.15 

0.39 

60 

0.06 

0.21 

0.00 

0.25 

0.23 

0.44 


80 

0.14 

0.31 

0.00 

0.36 

0.33 

0.59 


20 

0.01 

0.06 

0.00 

0.00 

0.00 

0.07 

SR2 

40 

0.03 

0.10 

0.00 

0.00 

0.00 

0.13 

60 

0.05 

0.12 

0.00 

0.00 

0.00 

0.17 


80 

0.06 

0.15 

0.00 

0.00 

0.00 

0.21 


20 

0.00 

0.41 

0.75 

0.00 

0.00 

0.75 

SR3 

40 

0.00 

0.51 

0.78 

0.00 

0.00 

0.78 

60 

0.00 

0.55 

0.86 

0.00 

0.00 

0.86 


80 

0.00 

0.64 

0.86 

0.00 

0.00 

0.86 


Table 1. The ratio of cross section that passes the trigger cut to the generated cross section for 200k 
generated events. Kinematic dependent tagging efficiencies are already incorporated into the cross 
sections. All leptons (e, //, r) are generated with a minimum pt > 10 GeV, but r decays to leptons 
can result in a p^^ <10 GeV. The columns in this table are not necessarily independent, as it is 
possible for an event to simultaneously pass multiple triggers. 


• pp ^ "i"" jZ + bb + ( 0 , l)j, 7*72^ ^ ; 

• pp^ -i*jz + 6(6) + (0,1,2)j, -i*jz 

• pp —>• 7*/Z + (0, l)2,3)j, 7*/Z —> ; 

• pp ^ + 66 + (0, l)j, 

• PP^W^ + 6(6) + (0,1,2)j, ^ 

• pp^w^ + (0,1,2,3)j, ^ 

• pp ^ W~^W~ + 66 + (0, l)j, ^ W~ —> 

• PP^ W+W- + 6(6) + (0,1, 2)j, W+ e+ue, W- 

• pp^W+W- + {0,1,2,3)j, W+ W- 

• PP^ ZW^ + 66 + (0, l)j, z ^ f+f-; 

• pp^ ZW^ + 6(6) + (0,1, 2)j, Z f+f"; 

• pp^ ZW^ + (0,1,2, 3)j, Z f+f-; 

where i = (e, /i, r) and j are light jets (ix, d, 5 , c, g) that can come from associated production. 
Each entry in the list above is produced with the number of quoted jets, and MLM matching 
and merging is incorporated to avoid double counting of the light jet production with the 
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Figure 4. Trigger ratios for each signal region, normalized to the produced and tagged cross section, 
based on varying the leading muon pt- The pCa trigger is based on a subleading electron pt = 12 
GeV, while /ie^ is based on a subleading electron pt = 17 GeV. The 2/i trigger for SR3 is based on 
a subleading muon pt = 15 GeV rather than the pt = 10 GeV discussed in the text, as the trigger 
response for the lower subleading pt is very similar to the single muon rate due to the minimum pt 
settings in the event generation stage and tagging thresholds. 


initial state radiation (MLM matching with XQCUT = 15 and QCUT = 20). The two 
largest contributions to our backgrounds are the inclusive Z production modes (included in 
the first three entries) and tt (included in the seventh entry), but these are effectively reduced 
by kinematic cuts. The kinematic distributions of the signal and backgrounds are included 
in Appendix A. 

Based on the kinematic distributions we examined, we have identified a number of possible 
cuts that improve the signal significance. These cuts are focused on reducing the ti and Z-\-nj 
backgrounds. The tt and other backgrounds with W^W~ lepton production can be reduced 
with cuts that involve the measurement, including a direct cut, as well as the transverse 
mass mp — 2p^^ ^ ^r(l ~ cos 9). Backgrounds with a Z resonance can be reduced by a 

cut on the dilepton invariant mass, In addition, a large fraction of the backgrounds 

producing both leptons and jets have a large total pp. Thus, we consider cuts on the scalar 
sum Hp — '"^^Pp + 

i b i 
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In the case of SR3, dilepton invariant mass cuts are implemented in a fixed range. While 
the branching ratio to dimuons is small (< 0.1%), the a branching 

ratio is similarly small, and the invariant mass peak of the direct decay is reconstructible with 
low smearing. Thus, it may be possible to observe the pseudoscalar with a resonance search 
methodology. For SR3, we consider only events within a 2-3 GeV invariant mass bin centred 
at the mass of the pseudoscalar. In contrast, the analysis for SRI and SR2 are based on a 
cut-and-count methodology, since the dilepton peak is significant smeared out due to the loss 
of information from the neutrinos originating from the r decays. For these signal regions, 
we do not employ a narrow invariant mass window and instead employ mu cuts to exclude 
backgrounds only. 

The cuts for SRI and SR2 are considered separately in each of two distinct scenarios: 
hard cuts are better for high luminosity searches and have a greater overall reach, while soft 
cuts are better for low luminosity search es. Kinematic threshold values fo r the considered 
cuts were chosen by maximizing (Jsig^L/^cisig * L + ai^kg * T + ^^sys^bkg * ^ systematic 

uncertainty of eggs — 0.2 and luminosity of L = 100/fb, while maintaining ~ 

0.5(0.8) for hard (soft) cuts for ma > 40 GeV. The dimuon signal region, SR3, is analyzed 
assuming only a single cut scenario, as background events with ^ ma generally have 
similar acceptance rates to the signal. 

The final cut values for each signal region are: 

• SRI hard: leading < 30 GeV, 12 < m^i < 35 GeV, Ht < 90 GeV, < 80 GeV. 

• SRI soft: leading < 40 GeV, 12 < ma < 45 GeV, Ht < 140 GeV, < 40 GeV, 

< 120 GeV. 

• SR2 hard: leading p^ < 40 GeV, 12 < m^^ < 45 GeV, Ht < 130 GeV, < 100 GeV. 

• SR2 soft: 12 < mu < 60 GeV, Ht < 190 GeV, mr < 45 GeV, < 140 GeV. 

• SR3: leading < 50 GeV, Ht < 120 GeV, < 120 GeV. 

The expected search reach using these cuts is given in the next section. Further details about 
the acceptance rates for each cut are provided in Appendix B. Alternative approaches for 
determining the cut regions, such as those incorporating repeated algorithmic refinements of 
the phase space, would optimize cuts for a single mass value and be unable to account for the 
full range of parameters we explore. Maximizing the acceptance rate for = 40 GeV would 
result in a poorer reach in gg values for ma — 80 GeV, for example. We feel our approach is 
more appropriate for a general search strategy. 

4 Results 

We can now investigate the extent to which the light pseudoscalar parameter space consistent 
with the Fermi signal can be probed by the searches we propose. 
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Figure 5. Discovery potential for the SRI signal region with hard (left) and soft (right) cuts for 
^sys = 0.1 (0.3) (top (bottom)). Hatched region is the region where no discovery is possible, regardless of 
luminosity, due to systematic uncertainties. The shading and labeled contours correspond to constant 
values of \og{L x fb) needed to achieve k = 3. 


Due to the low pseudoscalar mass region of interest in this study, as well as the cut-and- 
count search method for SRI and SR2, systematic uncertainties are a particularly challenging 
aspect of performing this search. To estimate the effect of systematic uncertainties, we con¬ 
sider two scenarios in addition to our two cut (hard/soft) scenarios - low systematics, with 
^sys — 10%, and high systematics, with Cgys = 30%. Our analysis of the discovery potential 
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Figure 6 . Discovery potential for the SR2 signal region with hard (left) and soft (right) cuts for 
^sys = 0.1 (0.3) (top (bottom)). Hatched region is the region where no discovery is possible, regardless of 
luminosity, due to systematic uncertainties. The shading and labeled contours correspond to constant 
values of \og{L x fb) needed to achieve k = 3. 


is based on a signal significance, given by 


k = 


Ns 

^Ns + Nh + e%,Nf 


(4.1) 


where Ng = ag * L and N), = (Tb * L are the number of signal and background events, 
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Figure 7. Discovery potential for the SR3 signal region for egys = 0.1 (0.3) (left (right)). Hatched 
region is the region where no discovery is possible, regardless of luminosity, due to systematic uncer¬ 
tainties. The shading and labeled contours correspond to constant values of log{L x fb) needed to 
achieve k = 3. 

respectively, after cuts for a given integrated luminosity, L. Contours of constant luminosity 
are plotted in Figures 5, 6 and 7. For small enough values of systematic uncertainties 
dominate the signal, and we expect that greater luminosity will be insufficient to illuminate 
any signal. Note that we have also verified that each signal data set considered has at least 
5 events after cuts. 

The soft cut scenarios of SRI and SR2 are optimized for early searches with low lumi¬ 
nosity, but suffer from a larger systematics-dominated region, since the total backgrounds 
are much larger. Thus, their ability to exclude the parameter space ends at approximately 
L = 10/fb integrated luminosity. Alternatively, hard cuts scenarios have a better reach with 
exclusions from L = 100/fb, though larger luminosity will be unlikely to push this boundary 
any further. 

As discussed, the expected sensitivites for each case are affected by three primary com¬ 
ponents: production, trigger and cuts. Production rates decrease with increasing mass, 
reducing the overall cross section and number of events at the LHC. In contrast to this, trig¬ 
ger response improves for heavier pseudoscalars, but has a significant effect on the lighter 
pseudoscalar scenario. However, the pseudoscalar is produced in association with b quarks, 
which results in a boost to the a that allows a large enough fraction of events to pass trigger 
and thereby make the search viable. Lastly, eliminating backgrounds resulting from the Z 
peak results in a choice of cut thresholds that has a larger impact on events from heavier 
pseudoscalar masses, especially for the hard cut scenarios. These issues combined result in 
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the typical shape observed in Figures 5 and 6, with reduced exclusion reach for both the 
lowest and highest mass scenarios. 

The dimuon search uses a different approach, incorporating a pseudo-resonance search 
methodology. While we do not fit a line-shape over the background and compare the signal, 
we employ a narrow invariant mass window with a sliding center that effectively estimates 
the result from such an approach. In practice, an approach that fits a line to the continuum 
background will reduce systematic uncertainties that are associated with the cut-and-count 
methodology, which requires simulations to estimate the backgrounds. As a result, we suspect 
that the low systematics scenario in Figure 7 is potentially a more realistic case, in contrast 
with the other signal regions, where low systematics may be overly optimistic. 

As a result of the relatively large width of the SM Z, combined with detector smearing 
effects, a dimuon resonance at 80 GeV will contend with increased backgrounds from the 
Z peak (which is why we do not consider heavier masses). If we assume similar systematic 
uncertainties for each signal type, then the most promising reach for the high rria region is 
in the lel/i signal regions, while the 1^1 r signal regions are more promising for the low rria 
regime. Note that the reach in the dimuon signal region is not as promising as the others for 
any part of the parameter space under the assumption of similar systematies. As mentioned, 
however, systematic uncertainties in the dimuon search will likely be smaller than in the other 
modes, and so all signal regions combine to form a complimentary and robust search strategy. 

Comparing Figures 1 and 5-7, we see that the searches we propose will cover a signif¬ 
icant portion of the otherwise unconstrained parameter space consistent with the Galactic 
Center excess in scenarios with light pseudoscalar mediators, even with rather low integrated 
luminosity. This region is both theoretically and phenomenologically well-motivated, and we 
encourage both ATLAS and CMS to consider searches along the lines of those presented here. 

5 Application to the NMSSM 

To illustrate the usefulness of our results in a UV-complete model, we can consider how our 
searches impact the Zs-symmetric NMSSM parameter space consistent with the excess. To 
set our conventions, we take the superpotential to be 

W = WMSSMl=o + XSHu-Hd+'^S^ (5.1) 

with soft supersymmetry breaking terms given by 

AVo/t = \Huf + \h£ + m| + XAxH^ ■ H^S + (5.2) 

Hatted quantities are chiral superfields. The lightest pseudoscalar mass eigenstate can be 
written in terms of the SU{2) and singlet pseudoscalar {A and a^, respectively) as 

a = A cos 9 -\- as sin 9 (5.3) 
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with effective couplings 

Qu — cosOcotp^ = cos0tan/3. (5.4) 

For sizable tan/3 and cos0 not too small, gd will be larger than 1. Our conventions follow 
those found in Refs. [23, 77], to which we refer the Reader for further details regarding the 
spectrum. 

There have been two scenarios proposed in the Zs-invariant NMSSM to explain the 
GC excess involving neutralino annihilation into SM particles through a light singlet-like 
pseudoscalar [16, 17] (see also Ref. [46] for an analysis of the general NMSSM, which in 
some cases may also be probed by the searches we present). The first involves a mixed 
singlino/Higgsino-like neutralino, which, to achieve a Standard Model- (and not singlet-) 
like 125 GeV Higgs, requires the lightest pseduoscalar to be a nearly pure singlet [16] (i.e. 
cos 0 <C 1). Since the pseudoscalar couplings to SM fermions are suppressed, to explain the GC 
excess this scenario requires 2m^ to within about a GeV precision, as well as additional 

Z-mediated contributions to the annihilation rate in the early universe to drive down the 
relic abundance. This would seem quite finely tuned, requiring a fortunate conspiracy of 
parameters to achieve. Instead, we focus on the second possibility, namely that the neutralino 
is bino/Higgsino-like. In this case, the singlet component of the 125 GeV Higgs is naturally 
small, and so the lightest pseudoscalar can feature a more significant amount of mixing 
between the singlet and SU{2) states. As a result, the requirement that the neutralino 
annihilation is on resonance is relaxed, allowing one to consider a much larger range of masses 
not precisely tuned to rria ~ 2m^ [16]. 

It is worth mentioning that analyses of the NMSSM subsequent to Ref. [16] have found 
somewhat different results, favoring the singlino/Higgsino scenario [17, 109]. However, taking 
the systematics into account in fitting the Fermi signal [12, 13], we find that the bino/Higgsino 
scenario is fully compatible with both the GC signal and the Fermi dwarf spheroidal limits. 
Another reason the bino/Higgsino scenario may have been disfavored in Ref. [109] is that the 
large pseudoscalar couplings to the SM fermions in the bino/Higgsino scenario are constrained 
by rare meson decays, in particular Bg ii~. As pointed out in Ref. [16], these constraints 
can be avoided rather straightforwardly by taking advantage of mild cancellations between 
the various SUSY contributions to BR{Bs /x+/i“). Such points can be difficult to sample 
in a large scan of the parameter space, as employed in Refs. [17, 109]. However, we have 
verified that the bino/Higgsino scenario is still in fact viable when taking these constraints 
into account, as claimed in Ref. [16]. 

The bino/Higgsino explanation for the GC excess maps directly onto our simplified model 
(only that the WIMP is a Majorana, instead of Dirac, fermion). To illustrate the effect of our 
searches on the viable bino/Higgsino parameter space of the NMSSM, we performed a Markov 
Chain Monte Carlo scan of the parameter space using NMSSMTools 4.4.0 [110], interfaced 
with micrOmegas 3.1 [111]. Motivated by the parameter space presented in Ref. [16], we 
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Figure 8. Application of our results to the Z 3 -symmetric NMSSM. The black (gray) contours corre¬ 
spond to the reach at 100 fb“^ (1 fb“^) for the hard (soft) cut scenarios and low systematics in the 
various search channels. The gray points are the result of a Markov Chain Monte Carlo scan of the 
parameter space (described in the text) consistent with all existing phenomenological constraints with 
no requirements on the LSP relic abundance or annihilation rate with parameters as in Eq. 5.5 and 
rriA = 550 GeV. The green, blue, and orange points correspond to points capable of explaining the 
Fermi signal and consistent with the recent dwarf spheroidal constraints for ttia = 500, 550, and 600 
GeV, respectively. The red band is an example of the NMSSM parameter space found to be consistent 
with the excess in Ref. [16]. The sample point of Table 2 below is indicated with a star. Note that it 
may be possible to choose parameters minimizing the haa coupling to fill in the rria < gd > ^ 

region, which we did not attempt in our scan. 


fixed 

A = 0.05, /i = 615 GeV, = 550 GeV, 

Ml = 45 GeV, M 2 = 1 TeV, M 3 = 2 TeV, (5.5) 

Mq 3 = M [/3 = 7.5 TeV, At = V6MQ^, 3 = 5-5 TeV 

with all other soft masses and triscalar couplings at 1 TeV, while varying tan/3, and A,^. 
We required all points to satisfy all existing constraints discussed earlier and implemented in 
NMSSMTools. The results of the scans are shown, along with our results for the LHC reach 
across the parameter space, in Fig. 8 . The gray points were generated without requiring 
the lightest supersymmetric particle (LSP) to explain the Galactic Center excess or satisfy 
constraints on its relic abundance. The green, blue, and orange points correspond to ttia = 
500, 550, 600 GeV and feature a bino-like LSP with a relic abundance compatible with WMAP 
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A 

K 


tan (3 

rriA 

9 

Ml 

M2 

0.05 

0.52 

-8.5 

21.8 

550 

615 

45 

1000 

run 

rria 


9d 

VLh? 

{av) [cm^/s] 

crsi [cm^] 

crsD [cm^] 

125.8 

67.5 

44.1 

16.4 

0.137 

1.48 X 10“26 

4.3 X 10“^® 

3.9 X 10“^^ 


Table 2. Example parameter space point in the NMSSM capable of explaining the GC excess and 
consistent with the Fermi dwarf spheroidal limits. All dimensionful parameters are in GeV unless 
otherwise stated. The remaining parameters are set to the values shown in Eq. 5.5. This point would 
likely be probed by the searches we propose at the 13 TeV LHC with 100 fb“^ of integrated luminosity. 


and Planck measurements (including a 2a theoretical uncertainty) [23] 

0.091 < < 0.138 (5.6) 


and compatible with both the Galactic Center excess and the dwarf constraints, 


1.0 < 


1 X 10“^® cm^/s 


< 1.5, 


(5.7) 


for self-conjugae dark matter. Points satisfying these constraints typically have small, but 
non-negligible, p-wave suppressed contributions at freeze-out, such as those involving the Z 
(but still consistent with limits on the invisible Z width). This slightly reduces the relic 
abundance relative to the value suggested by XX cl ^ bb annihilation alone and allows 
these points to circumvent the dSph limits. Note that we did not attempt to minimize the 
haa coupling, and so no points were found with 2ma < rrih and > 1. However, it might 
be possible to reach this parametric regime [23] as suggested in Ref. [16], whose results we 
show along with ours in Fig. 8 by the red band. These values were taken from Fig. 6 of 
Ref. [16] for = 35 GeV, while our scan was performed assuming ^ Mi — 45 GeV. 
Table 2 provides the detailed spectrum information for an example parameter space point 
consistent with the GC excess and which would be probed by a ^ at the 13 

TeV LHC. This point is marked by the black star in Fig. 8. Note also that our scan did not 
find points with > 18. Larger values of g^ are typically excluded by LHC limits on the 
heavy MSSM-like pseudoscalar for the values of tan [3 sampled. In theories that do not rely 
on mixing with the SM-like Higgs, these constraints, as well as those from h ^ aa decays, 
are often significantly relaxed or absent. 

The contours in Fig. 8 show the sensitivity of our proposed searches to the NMSSM 
parameter space consistent with the Galactic Center excess at both 1 fb“^ and 100 fb 
A significant portion of the favored region with sizable g^ would be probed by the 13 TeV 
LHC at these luminosities. Even more reach would be expected at the 14 TeV LHC. Our 
searches are complementary to h ^ aa observations as well as existing LHC searches for 
MSSM Higgs bosons and would access regions of the parameter space not currently probed 
by other experiments, providing a potential window into a dark sector difficult to access 
otherwise. 
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6 Summary and Conclusions 


Many dark matter models feature WIMPs that can be very difficult to observe at colliders. 
Scenarios of this type can be consistent with the Galactic Center excess observed by the 
Fermi Large Area Telescope. Exploring these “coy dark sectors” at the LHC suggests a shift 
away from missing transverse energy signals and towards direct signatures of the particle(s) 
mediating the interaction of the dark matter with the Standard Model. 

Models involving pseudoscalar mediators and consistent with the GC excess can be of 
the coy variety. A good fit to the Fermi signal can be provided by relatively light WIMPs 
annihilating through a pseudoscalar into b quarks. In many realistic scenarios this suggests 
substantial couplings of the mediator to down-type Standard Model fermions. The signal 
favors WIMP masses in excess of ^ 35 GeV, while current collider bounds often imply pseu¬ 
doscalar masses below 90 GeV (provided they satisfy constraints from LEP). An interesting 
and currently untested explanation of the GC signal thus involves a pseudoscalar with mass 
below about 90 GeV with sizable couplings to down-type fermions and small branching frac¬ 
tion into WIMPs. The latter is generically small in this scenario since the on-shell decay of 
the mediator into dark matter is often kinematically forbidden and because the pseudoscalar’s 
coupling to WIMPs is relatively small. Our study has attempted to extend LHC coverage to 
this scenario by taking advantage of the mediator’s enhanced couplings to Standard Model 
fermions (relative to those of a SM-like Higgs boson of the same mass) and studying the 
production and decays of the pseudoscalar involving down-type final states. 

To this end, we explored signals that include one to two 5-jets and with either r or /i 
lepton pairs in the final state. We employed a simplified model, in which we assumed that the 
couplings of the pseudoscalar to Standard Model fermions were proportional to their mass, 
modulo common scaling factors for down- and up-type fermions. While this need not be the 
case, this situation is common in UV completions involving Type H 2 Higgs doublet models, 
as in supersymmetry. Our results can be applied to models with different coupling structures 
by a straightforward re-scaling of the production cross-section and branching ratios. 

Due to the rather low pseudoscalar masses we consider, trigger is an important factor in 
the search reach. We thus performed an analysis of the trigger response of the signal, and 
explored cuts that were effective in improving the signal significance. Our search strategy 
comprises a signal excess analysis for the lelfi + 1 — 2b and lllr + l — 2b modes, including low 
luminosity (soft cuts) and high luminosity (hard cuts) scenarios, and a dilepton resonance 
search in the + 1 — 2b signal. Since signal excess searches suffer from large systematics 

from comparisons to simulated instead of data driven backgrounds, we also analyzed the 
impact of systematic uncertainties on the LHC reach in all three signal modes. 

In the most optimistic scenarios, we find that the LHC should be able to explore values 
of the reduced pseudoscalar coupling to down-type fermions as low as ^ 8 for 100/fb 
of integrated luminosity at ^/s =13 TeV. Even in more pessimistic scenarios with higher 
systematics, we find that the LHC should be able to explore down to gd ^ 10 for some 
values of m^. This reach, however, is highly dependent on the trigger settings, and so we 
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strongly recommend that the experimental collaborations attempt to account for this type 
of signal when finalizing their trigger thresholds for leptons, particularly those triggers for 
muons. The parameter space in the NMSSM not covered by /i ^ aa searches, with rria ^ 
60 — 80 GeV, should be explorable to some extent, and further optimization of the search 
strategy could focus on this narrow region of masses. More generally, the searches we propose 
are highly complementary to those already existing at the LHC or elsewhere, highlighting 
their importance in the interest of fully covering the parameter space in question. 

In summary, light pseudoscalars with significant couplings to Standard Model fermions 
are well-motivated mediators for dark matter annihilation and arise in many models, including 
those explaining the Fermi Galactic Center excess. In many cases, these new particles would 
have evaded previous searches but should be testable at the LHC. Significant regions of the 
parameter space can be explored even with low luminosity, and so this signal presents a 
possibility for ongoing examination throughout the full LHC program. 
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A Appendix: Kinematic Distributions 

The kinematic distributions of signals and backgrounds are included in this section. For 
the signals, 1.5 million events are generated for each of the three-body processes, a{pp 
bha)BR{a and two-body processes, a{pp b(h)a)BR{a For the back¬ 

grounds, each process discussed in Section 3 was generated with at least 1.5 million events, 
with some generated at higher multiplicity in order to achieve sufficient statistics to be confi¬ 
dent on the distributions. All distributions are generated from events remaining after applying 
trigger level cuts and then signal region tagging. Backgrounds are plotted additively, such 
that each successively larger background is added to the previous backgrounds. All distri¬ 
butions combine the generated events from multiple processes, such that the bin value for a 
single event is dependent on the specific process that resulted in the event generated. For 
this reason, two successive bins containing a single generated event each may have different 
weights. 
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Leading b-Jet (GeV) 




Figure 9. Leading 6-jet pr distribution. Trigger cuts and tagging are applied, but no other kinematic 
cuts are applied. From left to right, figures are for SRI (lel/i), SR2 {lilrh) and SR3 (2/i). 



Leading Lepton (GeV) 



Figure 10. Leading lepton (e, ja) pt distribution. For SR2, there is only one lepton. Trigger cuts 
and tagging are applied, but no other kinematic cuts are applied. From left to right, figures are for 
SRI (lel/i), SR2 llilrh) and SR3 (2/r). 





Figure 11. Sub-leading lepton pt distribution. For SR2, this is the r pr distribution. Trigger cuts 
and tagging are applied, but no other kinematic cuts are applied. From left to right, figures are for 
SRI (lel/i), SR2 llilTh) and SR3 (2/i). 
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Figure 12. Leading 6-jet pt distribution. Trigger cuts and tagging are applied, but no other kinematic 
cuts are applied. From left to right, figures are for SRI (lel/i), SR2 {lllTh) and SR3 (2/i). 




Figure 13. AR distribution between the leading 6-jet and leading lepton. Trigger cuts and tagging 
are applied, but no other kinematic cuts are applied. From left to right, figures are for SRI (lel/i), 
SR2 {lllTh) and SR3 (2/i). 




Figure 14. Dilepton invariant mass distribution. Trigger cuts and tagging are applied, but no other 
kinematic cuts are applied. For SR3 with /i+/i“ final states, the width of the a is 0(1) GeV, and thus 
the direct dimuon production peak is more pronounced than what is shown with 5 GeV bins. From 
left to right, figures are for SRI (lel/i), SR2 (lllrh) and SR3 (2/i). 
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Figure 15. Scalar sum of visible transverse momenta distribution. Trigger cuts and tagging are 
applied, but no other kinematic cuts are applied. From left to right, figures are for SRI (lel/i), SR2 
{lilTh) and SR3 (2/i). 



Figure 16. Transverse mass distribution. Trigger cuts and tagging are applied, but no other kinematic 
cuts are applied. From left to right, figures are for SRI (lel/i), SR2 {lilrh) and SR3 (2/i). 



50 100 150 200 250 300 

+?^(GeV) 



Figure 17. Scalar sum of lepton and missing transverse momenta distribution. Trigger cuts and 
tagging are applied, but no other kinematic cuts are applied. Prom left to right, figures are for SRI 
(lel/u), SR2 (UlTh) and SR3 (2/i). 
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B Cut Flow Matrices 


When examining the potential of enhancing the visibility of the signal through cuts, we 
considered a variety of possible kinematic variable distributions, some of which are shown in 
Appendix A. Of those examined, we chose to consider only those cuts in which the shape of 
the backgrounds was distinctly different than the shape of the signal for at least one of the 
signal regions so that cuts on the background had a larger fractional effect on the backgrounds 
than on the signal. 

The variables that most effectively improved the signal significance were Pt of the 
leading lepton (pt^), dilepton mass {mu), total scalar sum of visible momenta {Ht)^ transverse 
mass of the subleading lepton {mT 2 ndi)^ scalar sum of the lepton pT and {Wt)' 

Variables with components were most effective at eliminating backgrounds containing 
decays of W bosons, including mT 2 nM where backgrounds containing intermediate W’s have a 
longer tail on the distribution. Of note, we found that the transverse mass distribution based 
on the leading lepton pT had a longer tail for the signal, and so was not quite as effective. 

Since for example, is a component of multiple cuts, we consider the correlation 
between the events passing each pair of cuts in cut ffow matrices in Tables 3 through 7. 
Diagonal entries are the acceptance rate for the single cut labeled in both the column and 
row headers, where red text indicates background acceptance rates and black text indicates 
signal acceptance rates. Each off-diagonal entry in these tables represents the acceptance rate 
{A) of the cut labeled by the row (r) header on the events remaining after performing the 
cut in the column (c) header, such that each entry is given by 

Ar,c = Ar{Ac{(j))/Ac{(j). (B.l) 

For example, the upper right most entry of Table 3 shows an 87.4% acceptance rate for 
background events and 84.3% acceptance rate for signal events from applying the cut to 
the pool of events that already passed the^j. cut. The lower left most entry shows that 17.2% 
of background events and 87.7% of signal events pass the cut after applying the cut. 
Since the cut removes a similar number of events for the signal and background once the 

cut has been applied, the cut is superffuous once the cut has been applied, and 
thus should not be included in the final set of cuts. In fact, the distribution for signals 
and backgrounds have similar shapes once the cut has been applied, and thus no cut 
value will be effective. 
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PTe 

m££ 

Ht 

^T2ndt 

JflT 


0.130 

0.431 

0.158 

0.523 

0.115 

0.246 

0.335 

0.619 

0.272 

0.452 

0.874 

0.843 

PTe 

0.070 

0.596 

0.058 

0.490 

0.203 

0.483 

0.512 

0.742 

0.088 

0.477 

0.451 

0.689 

m££ 

0.076 

0.250 

0.305 

0.431 

0.087 

0.438 

0.204 

0.285 

0.084 

0.512 

0.185 

0.298 

Ht 

0.050 

0.537 

0.172 

0.565 

0.046 

0.243 

0.019 

0.373 

0.041 

0.336 

0.296 

0.631 

'^T2nM 

0.295 

0.582 

0.214 

0.539 

0.138 

0.648 

0.301 

0.499 

0.141 

0.554 

0.478 

0.538 

jflT 

0.172 

0.877 

0.200 

0.629 

0.055 

0.648 

0.390 

0.756 

0.087 

0.435 

0.026 

0.448 


Table 3. Cut flow matrix for SRI: lel/i+l —26+Oj signal with hard cuts. The cuts are: < 30 GeV, 

< 30 GeV, 12 < mu < 35 GeV, Ht < 90 GeV, < 20, < 80. The red/top entry in each 

cell is the acceptance rate for all backgrounds combined, while the black/bottom entry shows the 
acceptance rate for the signal. The total cross sections for the lel/i + 1 — 26 + Oj signal after applying 
the trigger cuts are ai^kg = 2187 fb and asig = 60.4 fb {rria = 60 GeV, ga = 25). 




PTe 

m££ 

Ht 

^T2nM 

T/ft 

Jflr 


0.312 

0.721 

0.324 

0.765 

0.310 

0.699 

0.400 

0.845 

0.497 

0.734 

0.809 

0.885 

PTe 

0.188 

0.886 

0.181 

0.835 

0.456 

0.851 

0.503 

0.943 

0.236 

0.833 

0.538 

0.920 

m££ 

0.161 

0.850 

0.408 

0.893 

0.162 

0.877 

0.320 

0.871 

0.187 

0.898 

0.336 

0.863 

Ht 

0.236 

0.909 

0.511 

0.876 

0.364 

0.771 

0.184 

0.776 

0.249 

0.764 

0.535 

0.923 

'^T2„m 

0.487 

0.912 

0.398 

0.893 

0.353 

0.918 

0.413 

0.883 

0.306 

0.896 

0.626 

0.897 

Tjri 

Jflr 

0.404 

0.970 

0.463 

0.870 

0.323 

0.918 

0.453 

0.941 

0.319 

0.791 

0.156 

0.790 


Table 4. Cut flow matrix for SRI: lel/i + 1 —26+Oj signal with soft cuts. The cuts are: < 50 GeV, 

< 40 GeV, 12 < m££ < 45 GeV, Ht < 140 GeV, < 40, < 120. The red/top entry in 

each cell is the acceptance rate for all backgrounds combined, while the black/bottom entry shows the 
acceptance rate for the signal. The total cross sections for the lel/r + 1 — 26 + Oj signal after applying 
the trigger cuts are cji^kg = 2187 fb and cisig = 60.4 fb {rria = 60 GeV, gd = 25). 
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PTe 

mu 

Ht 

^T2ndt 

JflT 


0.170 

0.396 

0.221 

0.441 

0.151 

0.267 

0.389 

0.488 

0.328 

0.320 

0.814 

0.715 

PTe 

0.339 

0.671 

0.262 

0.602 

0.517 

0.570 

0.644 

0.840 

0.403 

0.602 

0.723 

0.857 

mu 

0.063 

0.288 

0.140 

0.405 

0.071 

0.428 

0.086 

0.306 

0.087 

0.472 

0.104 

0.340 

Ht 

0.328 

0.596 

0.353 

0.673 

0.173 

0.345 

0.143 

0.483 

0.243 

0.408 

0.716 

0.808 

'^T2nM 

0.344 

0.401 

0.274 

0.495 

0.218 

0.546 

0.303 

0.419 

0.178 

0.496 

0.451 

0.395 

jflT 

0.373 

0.759 

0.215 

0.597 

0.114 

0.546 

0.389 

0.702 

0.197 

0.334 

0.078 

0.420 


Table 5. Cut flow matrix for SR2: 1-^lr+l —26+Oj signal with hard cuts. The cuts are: < 30 GeV, 

< 40 GeV, 12 < m« < 45 GeV, Hr < 130 GeV, < 25, < 100. The red/top entry in 

each cell is the acceptance rate for all backgrounds combined, while the black/bottom entry shows the 
acceptance rate for the signal. The total cross sections for the l^lr + 1 — 26 + Oj signal after applying 
the trigger cuts are cjijkg = 742 fb and Ggig = 84.3 fb {rria = 60 GeV, Qd = 25). 




PTe 

mu 

Ht 

^T2nM 

T/ft 

Jflr 


0.422 

0.778 

0.456 

0.821 

0.462 

0.762 

0.530 

0.874 

0.642 

0.776 

0.860 

0.912 

PTe 

0.564 

0.899 

0.521 

0.851 

0.844 

0.850 

0.717 

0.944 

0.693 

0.871 

0.817 

0.955 

mu 

0.198 

0.864 

0.293 

0.880 

0.181 

0.882 

0.257 

0.872 

0.255 

0.899 

0.334 

0.874 

Ht 

0.587 

0.878 

0.642 

0.867 

0.663 

0.772 

0.467 

0.781 

0.575 

0.795 

0.816 

0.957 

'^T2„m 

0.467 

0.766 

0.408 

0.786 

0.433 

0.782 

0.378 

0.781 

0.307 

0.767 

0.572 

0.777 

Tjri 

Jflr 

0.596 

0.889 

0.458 

0.851 

0.539 

0.782 

0.511 

0.930 

0.545 

0.768 

0.292 

0.759 


Table 6. Cut flow matrix for SR2: l^lr + 1 —26 + Oj signal with soft cuts. The cuts are: < 55 GeV, 

< 55 GeV, 12 < mu < 60 GeV, Ht < 190 GeV, < 45, < 140. The red/top entry in 

each cell is the acceptance rate for all backgrounds combined, while the black/bottom entry shows the 
acceptance rate for the signal. The total cross sections for the l^lr + 1 — 26 + Oj signal after applying 
the trigger cuts are a^kg = 742 fb and agig = 84.3 fb (m^ = 60 GeV, ga = 25). 
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PTe 

Ht 

^T2nM 

Tjrt 

J/Lt 


0.882 

0.888 

0.883 

0.902 

0.962 

0.950 

0.950 

0.885 

0.992 

0.962 

PTe 

0.317 

0.878 

0.317 

0.864 

0.770 

0.969 

0.349 

0.859 

0.546 

0.907 

Ht 

0.288 

0.815 

0.642 

0.855 

0.264 

0.762 

0.318 

0.767 

0.496 

0.853 

'^ T 2 „de 

0.688 

0.858 

0.703 

0.856 

0.767 

0.866 

0.639 

0.860 

0.818 

0.887 

jflT 

0.532 

0.947 

0.815 

0.918 

0.888 

0.979 

0.606 

0.901 

0.473 

0.874 


Table 7. Cut flow matrix for SR3; 2/u + 1 — 26 + Oj signal. The cuts are: < 60 GeV, < 50 GeV, 

Ht < 120 GeV, < 45, < 120. The red/top entry in each cell is the acceptance rate for all 

backgrounds combined, while the black/bottom entry shows the acceptance rate for the signal. The 
total cross sections for the 2/i + 1 — 26 + Oj signal after applying the trigger cuts are ai)kg = 7249 fb 
and asig = 108 fb {rria = 60 GeV, ga = 25). 
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C Variation of Exclusions 


As discussed in section 3, calculations of signal events using the 5FS is quite strongly de¬ 
pendent on the factorization and renormalization scales used. In MadGraphB, we employed a 
dynamic scale scheme that we then varied by an overall scaling factor between 0.5 and 2.0 
(see Eq. 3.1). This factor had the largest effect for low mass pseudoscalar calculations, with 
a factor of 0.5 reducing the total cross section by approximately 22% for = 20 GeV, while 
only reducing the total cross section by a factor of 4% at rria — 80 GeV. Alternatively, the 
authors of [103] use a fixed renormalization and factorization scale scheme based on the sum 
of the masses of the pseudoscalar and the on-shell b quark masses. Variations of this scale by 
a factor between 0.5 and 2.0 results in a cross section reduced by as much as 50%. 

In addition, our calculations of the backgrounds were performed at leading order. Higher 
order effects, as well as possible unaccounted-for experimental issues, may result in larger 
backgrounds than we predict. In order to address concerns regarding these two issues, we 
explore much more conservative contours determined by performing the same calculations 
but with a factor of 2.0 larger backgrounds, and separately with a factor of 0.5 smaller signal. 
Figures 18, 19 and 20 give these results. Of note, many regions of parameter space are still 
explorable at the LHC with 100/fb of integrated luminosity even in the more pessimistic 
scenarios. 
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Figure 18. Discovery potential contours for the SRI (lel/i) signal region with hard (left) and soft 
(right) cuts for Cgys = 0.1 (0.3) (top (bottom)) and conservative factors applied to the signal (dotted) 
and backgrounds (dashed) (solid lines show the original bounds without any factor applied to the 
signal or background). Contours correspond to constant values of log(L x fb) needed to achieve k = 3. 
The black lines represent the boundary of the systematics dominated region, the red lines represent the 
discovery potential at I/=10/fb, while the yellow lines represent the discovery potential for L=l/fb. 
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Figure 19. Discovery potential contours for the SR2 (l^lr) signal region with hard (left) and soft 
(right) cuts for esys = 0.1 (0.3) (top (bottom)) and conservative factors applied to the signal (dotted) 
and backgrounds (dashed) (solid lines show the original bounds without any factor applied to the 
signal or background). Contours correspond to constant values of log(L x fb) needed to achieve k = 3. 
The black lines represent the boundary of the systematics dominated region, the red lines represent the 
discovery potential at L =10/fb, while the yellow lines represent the discovery potential for L =l/fb. 
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Figure 20. Discovery potential contours for the SR3 (2/i) signal region for esys = 0.1 (0.3) (left (right)) 
and conservative factors applied to the signal (dotted) and backgrounds (dashed) (solid lines show the 
original bounds without any factor applied to the signal or background). Contours correspond to 
constant values of log{L x fb) needed to achieve k = 3. The black lines represent the boundary of the 
systematics dominated region, the red lines represent the discovery potential at L =10/fb, while the 
yellow lines represent the discovery potential for L =l/fb. 
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